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We present improved bounds on the CP-violating phase $ J J^ and 
on the decay-width difference Ar of the neutral B° s meson system ob- 
tained by the CDF experiment at the Tevatron collider . We use 6500 
B° — > J/ip4> decays collected by the dimuon trigger and reconstructed 
in a sample corresponding to integrated luminosity of 5.2 fb _1 . Besides 
exploiting a two-fold increase in statistics with respect to the previous 
measurement, several improvements have been introduced in the analy- 
sis including a fully data-driven flavour-tagging calibration and proper 
treatment of possible S-wave contributions. 
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1 Introduction 



Since the discovery of CP in 1964 in neutral kaon system pQ, CP violation plays 
crucial role in the development of standard model (SM) and probing a new physics 
(NP). In 1973 Kobayashi and Maskawa proposed, as one of the possible explanations 
for CP violation in kaon system, extension to six quarks model [2] in which the 
CP violation is explained through the quark mixing parametrised by the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix. A single complex phase in the CKM matrix is 
responsible for all CP violation in the SM. The observation of large CP violation in 
B° mesons by BABAR and Belle experiments [3] confirmed the SM and paved way to 
Nobel Prize award to Kobayashi and Maskawa in 2008. 

After confirmation of the SM focus shifted to search for a NP. One of the most 
promising processes is B s mixing governed by the CKM matrix element Vt s . The 
indirect information defines Vt s to be almost real, which translates to the fact that 
the CP violation due to the B s mixing is expected to be tiny in the SM. First 
measurements performed by the CDF and D0 experiments [I] showed about 1.5-2.0 
cr discrepancy with the SM, which caused large excitement in the community. In 
this proceedings we review most important aspects and results of the updated CDF 
analysis, using dataset corresponding to integrated luminosity of 5.2 fb _1 [5]. 

2 Candidate selection 

To select candidates we use an artificial neural network (ANN) trained on a signal 
sample from simulation and data events from a B s mass sidebands as background. 
Inputs to the ANN are transverse momenta of B s and <fi mesons, particle identifica- 
tion for kaons and muons and quality of the kinematical fits to the candidates. For 
the result presented here we choose requirement on the ANN output which minimises 
uncertainties on the CP violating phase f3j^. The best point is found by performing 
simulated experiments with three different true values of fij^ (0.02, 0.3 and 0.5) 
and single value for Ar and amplitudes with amount of signal and background corre- 
sponding to different requirements on the ANN output and selecting the value which 
provides smallest parabolic uncertainty on the fi J J^$. In Fig. [I] we show an example 
of such simulation for true $ J J^ value at the SM expectation. As result we select 
candidates with the ANN output larger than 0.2, which is less stringent requirement 
than one used in previously. The resulting invariant mass distribution contains about 
6500 signal events and is shown in Fig. [TJ 

3 Flavour tagging 

The flavour tagging is one of the important tools in the analysis. Its task is to deter- 
mine whether reconstructed candidate was produced as B s or B s . Two algorithms are 
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Figure 1: An example of the dependence of parabolic uncertainty of f3 J J^^ on ANN 
output requirement (left). The invariant mass distribution of the selected candidates 
used in the analysis (right). 



used at hadron colliders. The first one, called opposite side tagging, explores the fact 
that 6-quarks are dominantly produced in bb pairs, so determination of the flavour 
of other 6-hadron determines also the flavour of signal one. The second algorithm, 
called same side tagging, explores properties of hadrons produced in hadronization of 
b quark into B s mesons. 

The opposite side algorithm explores several sources of information from the non- 
signal 6-hadron in an event. The most clean explores the fact that about 10% of 
6-hadrons decay to a final state containing lepton. At CDF we use electrons and 
muons and given that most of them arise from the 6-hadron decays in events with 
reconstructed B s , they provide clean information, but at expanse of small efficiency. 
The second source uses fact the most abundant b — > c — > s quark level transition 
yielding a final state which often contains charged kaon, which determines the flavour 
of the 6-hadron. The most efficient source identifies jet containing b quark and cal- 
culates weighted charge of the tracks in jet to determine, whether original quark had 
positive or negative charge and thus determine flavour. While this algorithm is ef- 
ficient, its purity is rather small. We combine all three sources into single decision 
using neural network. The quality of decision is predicted for each event and checked 
on the fully reconstructed B + events. The overall performance of the algorithm is 
about 1.2%, which can be understood as having 1.2% of the overall statistics with 
perfectly known production flavour. 

The same side algorithm exploits fragmentation process. To form a B s out of b 
quark we need to attach s quark to it. The s quark comes from a pair generated out 
of the vacuum and after forming B s an s quark remains to form other hadron. If it 
ends up in the charged kaon it can be used to determine the production flavour. The 
challenge is to identify right kaon in an event with many tracks. Another challenging 
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Figure 2: The invariant mass distribution of B s — > D s n with D s — > <f)ir (left) and the 
amplitude scan in the full dataset (right). 



part is to calibrate quality of the decisions on a data as only process available is 
B s mixing. Despite all challenges, CDF uses this algorithm with good success. In 
calibration we use data spanning same period as our dimuon dataset. For calibration 
we reconstruct about 12900 B s — > D s (3)tt signal events. The predicted quality of 
decision for each event is scaled by a global scaling factor which is determined in the 
fit for B s oscillations. We show the invariant mass distribution of channel contributing 
about half of the statistics together with the amplitude scan on full sample in Fig. |2j 
The resulting tagging performance is about 3.2%. The obtained B s mixing frequency 
is Am s = 17.79 ±0.07 ps _1 with uncertainty being statistical only. Full details of this 
calibration can be found in Ref. [5]. 



4 Fit description 

On the limited space available we cannot describe all details of the fit, so we will touch 
only main features with some emphasis on changes compared to the previous versions 
of the analysis. The full details are spelled out in Ref. [6]. As spin zero B s decays 
into two spin one particles {J/ip and 4>) three different amplitudes corresponding to 
different angular momentum are involved. In CDF we use a basis in which three 
amplitudes are written in terms of polarisation amplitudes. The three amplitudes 
give six angular terms, three terms corresponding to the squares of amplitudes and 
three corresponding to the interference between amplitudes. Each of the six terms 
has its own angular and decay time dependence. Some of the terms exhibit usual 
sin(Am s t) behaviour and some are proportional to cos(Am s t). Thanks to the non- 
zero decay width difference (Ar) in the B s system, depending on the size of the 
Ar and polarisation amplitudes one can gain considerable information on the CP 
violating phase also without resolving the oscillation pattern. 
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Figure 3: The contours in /3/^-Ar plane after full coverage adjustment (left). The 
likelihood profile for the fraction of s-wave contribution (right). 

Since the first analysis there is an ongoing discussion whether there can be con- 
tribution from non- resonant J/ipK + K~ or J/tp /o (980) with /o(980) — > K + K~ decays, 
collectively named s-wave. Those were neglected in previously and there is an open 
question about the size of possible contribution and whether it can bias result. In the 
presented result we introduce an s-wave component which is allowed to float during 
the analysis. This additional component introduces four new terms, one correspond- 
ing to the square of the s-wave amplitude and three for the interference between 
s-wave contribution and decay to J/ip<p. 

5 Results 

Similar to the previous iterations, the likelihood in case of all parameters floating 
still show non-Gaussian behaviour. In order to deal with the likelihood behaviour 
and ensure well defined statistical meaning of the result we construct confidence 
regions m the plane, which are shown in Fig. [3j For the standard model 

value of and Ar we obtain the p-value of 44% which corresponds to about 0.8 
standard deviations. Minimising also over Ar we obtain the p-value for standard 
model of 31% with #W e [0.02,0.52] U [1.08,1.55] at 68% confidence level and 
pW* e [-tt/2, -1.44] U [-0.13, 0.68] U [0.89, ir/2) at 95% confidence level. The amount 
of s-wave is consistent with zero and the likelihood profile over parameter describing 
its amount is shown in Fig. |3j 

In addition we also perform a fit in which we fix fi J J^ = 0. Such configuration 
provides well behaving likelihood allowing to provide point estimates of several inter- 
esting quantities even if the result should be interpreted in the context of standard 
model. From this fit we measure 

CT = 458.6 ± 7.6(stat) ± 3.6(syst) /im, 
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Ar = 0.075 ± 0.035(stat) ± O.Ol(syst) ps~\ 

|A||| 2 = 0.231 ±0.014(stat) ± 0.015(syst), 

\A \ 2 = 0.524 ± 0.013(stat) ± 0.015(syst), 

<f> ± = 2.95 ± 0.64(stat) ±0.07(syst). 

In all cases, those are the most precise measurements of those quantities up to date. 
The strong phase 6\\ is close to the symmetry point at ir which makes estimate of its 

value and uncertainty unreliable and therefore we do not provide result for it. 

6 Conclusions 

We presented the updated measurement of the CP violating phase fi J J^ in B s — > J/tp(f) 
decays from CDF experiment. Using 5.2 fb _1 of data we obtain bounds which are 
significantly stronger than our previous results. The data itself are consistent with 
the standard model at 0.8 standard deviation level. The improvement compared to 
the previous result is better than just simple scaling by the amount of collected data. 
Few possible improvements are still available in addition to collecting more data. On 
the data size itself, we expect to double our dataset by the end of 2011 with ongoing 
discussion for another 3 years extension to Tevatron running. 
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